The nucleon electromagnetic form factors from Lattice QCD 
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We evaluate the isovector nucleon electromagnetic form factors in quenched and unquenched QCD 
on the lattice using Wilson fermions. In the quenched theory we use a lattice of spatial size 3 fm at 
P = 6.0 enabling us to reach low momentum transfers and a lowest pion mass of about 400 MeV. 
In the unquenched theory we use two degenerate flavors of dynamical Wilson fermions on a lattice 
of spatial size 1.9 fm at /3 = 5.6 and lowest pion mass of about 380 MeV enabling comparison with 
the results obtained in the quenched theory. We find that unquenching effects are small for the 
pion masses considered in this work. We compare our lattice results to the isovector part of the 
experimentally measured form factors. 

PACS numbers: 11.15.Ha, 12.38.Gc, 12.38.Aw, 12.38.-t, 14.70.Dj 



I. INTRODUCTION 

The elastic nucleon electromagnetic form factors are 
fundamental quantities characterizing important features 
of neutron and proton structure that include their size, 
charge distribution and magnetization. An accurate de- 
termination of these quantities in lattice QCD is timely 
and important because of a new generation of precise 
experiments. In particular, polarization experiments |lj 
that measure directly the ratio of the proton electric to 
magnetic form factor, [i p G E jG v M , have shown a qualita- 
tive different behavior than the traditional Rosenbluth 
separation. The ratio fi p G p E /G p M instead of being ap- 
proximately constant falls off almost linearly with the 
momentum transfer squared, q 2 . This means that the 
electric form factor falls faster than the magnetic. Re- 
cent reviews on the experimental situation can be found 
m Refs. Hi. Precise lattice data for the nucleon form 
factors for large values of q 2 will enable comparison with 
experiment and could lead to an understanding of the 
approach to asymptotic scaling. Furthermore, access to 
low momentum transfers will enable a better determina- 
tion of phenomenologically interesting quantities such as 
the root mean squared (r.m.s.) radius of the transverse 
quark distribution in the nucleon 0. To access small 
momentum transfers we need a lattice with large spa- 
tial extent, L, since the smallest available momentum is 
2-k/L. Although large momentum transfers are in princi- 
ple available on typical lattices, the Fourier transform of 
two and three point functions becomes noise-dominated 
for momentum transfers beyond about 2 GeV 2 , limiting 
the range of high q 2 values that can be extracted accu- 
rately. 

In this work we calculate the isovector nucleon form 
factors as a function of the momentum transfer in lattice 
QCD both in the quenched approximation and with two 
dynamical Wilson fermions. A recent study of the nu- 
cleon form factors was carried out in the quenched theory 
using improved Wilson fermions The current work 



builds on the ingredients of the previous lattice calcu- 
lation and obtains results with higher accuracy at lower 
momentum transfers and pion masses. This enables us to 
determine the momentum dependence of the form factors 
accurately enough to have a meaningful comparison with 
experiment. A number of phenomenologically interest- 
ing quantities such as the r.m.s radii and dipole masses 
are extracted. Furthermore, we improve the accuracy 
of the results by constructing an optimal source for the 
nucleon allowing the maximum number of lattice mo- 
mentum vectors to contribute. The two form factors are 
then extracted using an overconstrained analysis that in- 
cludes all possible lattice measurements for a given value 
of q 2 . For the quenched calculation we use a lattice of 
size 32 3 x 64 at (3 = 6.0, which corresponds to a lattice 
spacing a ~ 0.09 fm, obtained either by using the nucleon 
mass at the chiral limit or the Sommer scale. In order 
to assess quenching effects, we also evaluate these form 
factors in the unquenched theory using dynamical Wil- 
son configurations simulated for quark masses that give 
pions of mass 690 MeV and 509 MeV on a lattice of size 
24 3 x 40 @ and 380 MeV on a lattice of size 24 3 x 32 at 
(3 = 5.6. The lattice spacing is about 0.08 fm determined 
from the nucleon mass at the chiral limit. This value 
is consistent with the value extracted using the Sommer 
scale over the range of quark masses used in this work. 

In lattice QCD, elastic matrix elements involving 
one-photon exchange require the evaluation of three- 
point functions. The standard procedure to evaluate 
three-point functions like the one we need here, namely 
G N i N (t%,ti; p ', p; T), is to compute the sequential 
propagator. This can be done in two ways: In an early pi- 
oneering work, where matrix elements of a number of dif- 
ferent hadronic states were evaluated @, the method of 
choice was to couple the photon to a quark at a fixed time 
t\ carrying a fixed momentum q. Within this scheme 
the form factors can only be evaluated at one value of 
the momentum transfer. Since the current must have a 
fixed direction and a fixed momentum, this approach is 
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referred to as the fixed current approach. This method 
allows one to use any initial and final state without re- 
quiring further inversions, which are the time consum- 
ing part of the evaluation of three-point functions. In 
the second approach, which is the method used in recent 
studies [MEi IH1> one requires that the initial state, cre- 
ated at time zero, and the final state, annihilated at a 
later fixed time t 2 , have the nucleon quantum numbers. 
The current can couple to any intermediate time slice t± 
carrying any possible value of the lattice momentum and 
having any direction. Therefore, within this scheme, with 
a single sequential propagator, one is able to evaluate all 
possible momentum transfers and current orientations. 
Since the quantum numbers of the final state are fixed, 
we refer to the second method as the fixed sink method. 
Clearly the fixed sink method is superior if our goal is the 
accurate determination of the momentum dependence of 
the nucleon form factors. 



II. LATTICE TECHNIQUES 

The nucleon electromagnetic matrix element for real 
or virtual photons can be written in the form 



N(p',s') \j„\N(p,s)) 



M 2 N 



E N (p>) E N (p) 



1/2 



u(p', s')Ofj,u{p, s) , (1) 



where p(s) and p'(s') denote initial and final momenta 
(spins) and Mn is the nucleon mass. The operator 
can be decomposed in terms of the Dirac form factors as 



2M 



N 



(2) 



where -Fi(O) = 1 for the proton since we have a con- 
served current and F 2 (0) measures the anomalous mag- 
netic moment. They are connected to the electric, Ge, 
and magnetic, Gm, Sachs form factors by the relations 



G M (q 2 ) = F^q 2 ) + F 2 (q 2 ) 



F 2 {q 2 ) 



(3) 



To extract the nucleon matrix element from lattice mea- 
surements, we calculate, besides the three point function 
G i (£2, ti\ p ', p; r), the nucleon two-point function, 
G NN (t,p), and look for a plateau in the large Euclidean 
time behavior of the ratio 



H(t 2 ,ti;p',p 



<G^(t 2 , p';r 4 )> 



(G^(t 2 -t 1 ,p;r 4 )) (G NN (t 1 ,p';T 4 )) (G NN (t 2 ,p';T 4 )) 
(G^(t 2 -tx,p';r 4 )} (G^(ti,p;r 4 )> <G™(t 2 , P ; r 4 )> 

=>■ n(p ,p;T;n) . 
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(4) 



We use the lattice conserved electromagnetic current, 
j^{x), symmetrized on site x by taking 

f( x )^\j»( x )+j»(x-{j)]/2 (5) 

and projection matrices for the Dirac indices 
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Throughout this work we use kinematics where the fi- 
nal nucleon state is produced at rest and therefore q = 
p' p = —p. Since we aim at obtaining the full q 2 
dependence of the form factors, we evaluate the three 
point functions with sequential inversions through the 
sink. We fix t 2 = 11(12) in lattice units for the quenched 
(unquenched) Wilson lattices and search for a plateau 
of R(t 2 ,t 1 - 1 p' ,p ;T;n) as a function of t\. Q 2 = —q 2 
denotes the Euclidean momentum transfer squared. 

We can extract the two Sachs form factors from the ra- 
tio of Eq. (gj by choosing appropriate combinations of the 
direction fi of the electromagnetic current and projection 
matrices T. Provided the Euclidean times t\ and t 2 — t\ 
are large enough to filter the nucleon ground state, the 
ratio becomes time independent. Inclusion of hadronic 
states in the two- and three-point functions leads to the 
expressions written in Euclidean space 



n(o,- q; r fc \n = %) = c 



2M N 



n(0,-q;r 4 ;/* = i) = c 



jk qj Gm(Q 2 



G e (Q 2 ) 



2M 



N 



n(o,- q; r 4 ;/ i = 4) 



c En + Mn_ Ge{q2) 



2M 



N 



(7) 



(8) 



(9) 



where C 



2Mi 



is a factor due to the normal- 



-E'jvC-E'jv+Miv) 

ization of the lattice states. The first observation regard- 
ing these expressions is that the polarized matrix element 
given in Eq. l(7)l. from which the magnetic form factor is 
determined, does not contribute for all momenta q. In 
the lattice study of the 7 TV — > A transition [lCl , we 
dealt with a similar situation where the naive A source 
was not optimal in the sense that, with one sequential 
propagator, not all lattice momentum vectors resulting in 
the same value of Q 2 contributed and an optimal source 
for the A was needed. Similarly, here one can construct 
an optimal linear combination for the nucleon sink that 
leads to 

3 C { 

S m (q]i) = ^II(-q ;T k ;fi = i) = ^rr-i (>2 - Pz)hA 

k=i N 
+(ps -Px)h,i + (pi-P2)5 3 , 1 }g m {Q 2 )(10) 
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and provides the maximal set of lattice measurements 
from which Gm can be extracted requiring one sequen- 
tial inversion. No such improvement is necessary for the 
unpolarized matrix elements given in Eqs. (HJ and ©, 
which yield Ge with an additional sequential inversion. 

Unlike the 7 TV — > A transition, the 7 N — > N transi- 
tion contains isoscalar photon contributions. This means 
that disconnected loop diagrams also contribute. These 
are generally difficult to evaluate accurately since the all- 
to-all quark propagator is required. In order to avoid 
disconnected diagrams, we calculate the isovector form 
factors. Assuming SU(2) isospin symmetry, it follows 
that 



2 1 - 



2 1 - 

( n\(-u^u - -d-y^d^n) 

(p\{uY'U-dj fi d)\p). (11) 



One can therefore calculate directly the three-point func- 
tion related to the right hand side of the above relation 
which provides the isovector nucleoli form factors 



G E (q 2 ) = G P E (q 2 ) — G E (q 2 ), 
GM(g 2 ) = GiM^-GlM 2 )- 



(12) 



The isovector electric form factor, Ge, can therefore be 
obtained from the connected diagram considering either 
the spatial components of the electromagnetic current 
as given in Eq. (JSJ or the temporal component given in 
Eq. 10, while Eq. Ill Oil is used for the extraction of the 
isovector magnetic form factor, Gm- 

Besides using an optimal nucleon source, the other im- 
portant ingredient in the extraction of the form factors is 
to take into account simultaneously in our analysis all the 
lattice momentum vectors that contribute to a given Q 2 . 
This is done by solving the overcomplete set of equations 



P(q; A1 )=73(q;/i).F(Q 2 ) 



(13) 



the different sink types, F — 



and D is an M x 



where -P(q; /i) are the lattice measurements of the ratio 
given in Eq. Q having statistical errors and using 

G E 
Gm 

2 matrix which depends on kinematical factors with M 
being the number of current directions and momentum 
vectors contributing to a given Q 2 . We extract the form 
factors by minimizing 



N 

El 

fc=i 



tl'k 



(14) 



using the singular value decomposition of D. Given the 
fact that one can have a few hundred lattice momentum 
vectors contributing in the evaluation of the form factors, 
the statistical precision is highly improved. Phenomeno- 
logically interesting quantities like the r.m.s. radii and 
magnetic moments can thus be obtained with increased 
precision. There is an additional advantage by includ- 
ing momentum transfers q as well as — q in our analysis. 
The lattice conserved current given in Eq. JSJ differs from 



the local electromagnetic current ^{x)^^ip{x) by terms 
of 0(a). However when we average over q and — q these 
0(a) terms vanish. 



E 
o 



local— local 
- x local— smeared 



.5 - 



a local— smeared with h 
* local— smeared with A! 




x smeared— smeared 
A smeared— smeared with hyp 
* smeared— smeared with APE 



10 



15 

t/a 



20 



25 



30 



FIG. 1: The rho effective mass as a function of the time 
separation on a 16 3 x 32 quenched lattice at f3 = 6.0 and 
ft = 0.153 using Dirichlet boundary conditions in the tempo- 
ral direction. In the upper graph, filled triangles show results 
obtained with local source and sink, crosses with Wupper- 
tal smeared source and local sink and open triangles (aster- 
isks) with Wuppertal smeared source using hypercubic (APE) 
smearing for the gauge links used in the construction of the 
hopping matrix i/(x, z; U(t)). The lower graph shows with 
crosses results obtained using Wuppertal smeared source and 
sink and with open triangles (asterisks) results with Wupper- 
tal smeared source and sink where hypercubic (APE) smear- 
ing is applied to the spatial links entering the hopping matrix. 



Smearing techniques are routinely used for achiev- 
ing ground state dominance before the signal from the 
time correlators is lost in the noisy large time limit. 
We use gauge invariant Wuppertal smearing, d{x, t) — > 
d smcar (x, t), at the source and the sink. We smear the 
fermion interpolating fields according to ^2 



(x, t) = F(x, z; U(t))d(z, t) 



(15) 



with the gauge invariant smearing function constructed 
from the hopping matrix H: 



F(x,z;U(t)) = (l + aH) n (x,z;U(t)), 



(16) 
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pertal smearing both to the source and sink, we see a 
reduction in the gauge noise in particular when hyper- 
cubic smearing is used. There is little effect on improv- 
ing ground state dominance since Wuppertal smearing 
on both source and sink very effectively cuts down ex- 
cited state contributions already after a time separation 
of a couple of time slices. Given the better noise re- 
duction observed when we use hypercubic smearing we 
choose to apply this smearing to the links that enter the 
hopping matrix H. For the parameters that enter the 
hypercubic smearing we use the same ones as those of 
Ref. [l^ . The parameters for the Wuppertal smearing are 
then optimized so that ground state dominance for the 
nucleon is optimal. We find that the parameters a — 4 
and n — 50 produce optimal results. These values are 
the same as those obtained without applying hypercu- 
bic smearing. Whereas Wuppertal smearing is applied to 
the source and the sink in all our computations to ensure 
ground state dominance at the time slice of the insertion 
of the electromagnetic current, hypercubic smearing is 
only done in the case of the unquenchcd configurations. 
This is because self averaging is less effective on smaller 
lattices causing the gauge noise to be more severe in the 
unquenched case where the simulations were done on a 
smaller lattice. 



FIG. 2: The nucleon effective mass as a function of the time 
separation on a 16 3 x 32 quenched lattice at (3 = 6.0 and 
ft = 0.153. The notation is the same as that of Fig. Q 



where 



£T(x, z; U(t)) = J2 ( *)<W-i + Uj{x - i, i)<5 x , 



y+i 
. 

It is well known that smearing introduces gauge noise in- 
creasing the errors on the extracted effective masses in 
particular when Wuppertal smearing is applied to both 
source and sink. An efficient way to reduce the the ul- 
traviolet fluctuations is to smooth the gauge fields at 
the time slice of the source or the sink where Wuppertal 
smearing is carried out. One can apply various smooth- 
ing techniques such as APE 13|, stout Q or hypercu- 
bic smearing on the gauge fields that are used in the 
hopping matrix. We found that both APE and hyper- 
cubic smearing reduce the noise and at the same time 
improve further the ground state overlap. In Figs. ^ and 
121 we show the effective mass for the rho meson and the 
nucleon respectively, using a lattice of size 16 3 x 32 with 
Dirichlet boundary conditions in the temporal direction 
to utilize the full time extent of the lattice. When only 
the source is smeared, both APE and hypercubic smear- 
ing improve the ground state overlap to such an extent 
that the plateau value is reached within a time separa- 
tion as short as two time slices. When we apply Wup- 



III. DISCUSSION OF EXPERIMENTAL 
RESULTS 

As explained in section II, we only compute the isovec- 
tor part of the form factors given in Eq. (|12|l . Therefore, 
to compare with experiment, it is necessary to extract 
from the experimentally available proton and neutron 
data the isovector contribution. In order to do this we 
need to interpolate the proton and neutron data to the 
same Q 2 values. 
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FIG. 3: The isovector electric form factor, Ge, extracted by 
interpolation from the measured proton and neutron electric 
form factors. 



In Fig. [3] we show the 
electric form factor 



tne proton and neutron data for the 
[lrl firj l . As can be seen, we need to 
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FIG. 4: The isovector magnetic form factor, Gm, extracted by 
interpolation from the measured proton and neutron magnetic 
form factors. 
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FIG. 5: The ratio of isovector form factors Ge over Gm as 
compared to the corresponding ratio of proton form factors 
from recent polarization experiments 



extrapolate the neutron electric form factor G E at low 
momentum transfers and the proton electric form factor 
G P E at intermediate Q 2 values in the range 0.25 < Q 2 < 
1.0 GeV. In order to interpolate the neutron data we 
consider the Galster parametrization 0] 



G N E {Q 2 ) = 



1 + 5.6r 



G d (Q 2 ) , 



(18) 



where t = Q 2 /AM N 2 , G d {q 2 ) = 1/(1 + Q 2 /0.71) 2 and 
fj, n = —1.91315, which provides a good description of 
the data. We calculate the derivative needed for the 
interpolation between measured data using the Galster 
parametrization. Similarly, in order to interpolate the 
proton data, we fit to a dipole form and then use the 
fitted form to compute the derivative needed for the in- 
terpolation. Having G E and G P E at the same value of 
Q 2 , we then find the isovector contribution via Eq. H12|) 
and plot the resulting isovector Ge in Fig. |21 As ex- 
pected, the difference between G P E and the isovector part 
is small due to the smallness of G E . A similar analysis 



is done for the magnetic form factor using the data of 
Ref. 0, 0| . The derivative needed for the interpola- 
tion is computed either from the best dipole fit function 
to the available data or, in the cases where we have two 
measurements close to the value of Q 2 that we are inter- 
ested in, by using a finite difference approximation to the 
derivative. The resulting isovector magnetic form factor 
is shown in Fig. 0] In Fig. [3] we compare the extracted 
isovector ratio [iGe/Gm to recent measurements of the 
proton ratio U wG p F / 'G p , that showed an unexpected Q 2 
dependence 0, Il9j |. The isovector ratio for Q 2 < 1 GeV 2 
decreases faster with Q 2 than /i p G p E /G p M , whereas for 
Q 2 > 1 GeV 2 it remains approximately constant. One of 
the goals is to compare this behavior with lattice calcu- 
lations. 



IV. LATTICE RESULTS 

As pointed out in the Introduction, the purpose of this 
work is to obtain accurate results over a large range of 
momentum transfers. For this reason the quenched calcu- 
lation is done on a lattice of size 32 3 x 64 enabling us to 
reach momentum transfers as low as about 0.15 GeV 2 . 
The highest momentum transfer that is accessible at 
(3 = 6.0 is 27r/a ~ 13 GeV. However statistical errors do 
not allow us to reach this maximum value. For the un- 
quenched calculation we use configurations generated by 
the SESAM collaboration |ona lattice of size 24 3 x 40 
and the DESY-Zeuthen group Q on a lattice of size 
24 3 x32at/3 = 5.6. At this value of (3 the lattice spacing 
is close enough to the lattice spacing of the quenched lat- 
tice so that finite a-effects are comparable. Differences 
between the two evaluations can then be attributed to 
unquenching effects. In Table [I] we give the parameters 
of our calculation. 

The lattice spacing is determined from the mass of 
the nucleon in the chiral limit. We use two different 
Ansdtze for extrapolating to the chiral limit: one is 
aM N = aM N (0)+c o ml and the other aM N = aM N (0) + 
cim 2 , + C2"m\. This provides an estimate for the system- 
atic error in the extrapolated value of the nucleon mass 
which dominates the overall error quoted in Table [I] for 
the lattice spacing. We would like to point out that using 
the Sommer scale to set the scale and taking tq = 0.5 fm, 
we find values that are consistent with the ones extracted 
from the nucleon mass. For the quenched lattice us- 
ing for ro/a the values given in Ref. |20j| . we obtain 
a- 1 = 2.15 GeV (a = 0.093 fm). For the unquenched lat- 
tice, the same definition gives at K, = 0.1575, k = 0.1580 
and k = 0.15825 a" 1 = 2.42(4) GeV, a" 1 = 2.47(3) GeV 
and a -1 = 2.56(6) GeV respectively. These values 
are consistent with the value of a -1 = 2.56(10) GeV ex- 
tracted from the nucleon mass at the chiral limit. The 
choice of the lattice spacing affects the physical value of 
Q 2 but not the values of Ge and Gm since they are com- 
puted in dimensionless units. In other words, a change 
in the lattice spacing a stretches the curves for Ge and 
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TABLE I: In the first column we give the number of config- 
urations and in the second column the value of the hopping 
parameter, re, that fixes the bare quark mass. In the third and 
fourth columns we give the pion and nucleon mass in lattice 
units. The values of the lattice spacing a are determined from 
the mass of the nucleon at the chiral limit. The unquenched 
configurations at re = 0.1575 and re = 0.1580 are provided 
by the SESAM collaboration and at k = 0.15825 by the 
DESY-Zeuthen group Q. 



number of confs 


re 


am-n 


clMn 


Quenched 32 3 


x 64 a~ 


1 = 2.14(6) GeV 


200 


0.1554 


0.263(2) 


0.592(5) 


200 


0.1558 


0.229(2) 


0.556(6) 


200 


0.1562 


0.192(2) 


0.518(6) 




re c =0.1571 


0. 


0.439(4) 


Unquenched 24' 


! x 40 a 


1 = 2.56(10) GeV 


185 


0.1575 


0.270(3) 


0.580(7) 


157 


0.1580 


0.199(3) 


0.500(10) 


Unquenched 24' 


' x 32 a 


1 = 2.56(10) GeV 


200 


0.15825 


0.150(3) 


0.423(7) 




re c = 0.1585 


0. 


0.366(13) 



Gm along the Q 2 axis and thus changes the slope. 

All the quenched results for the form factors are ob- 
tained using 200 configurations and three values of the 
hopping parameter re. The values of re chosen are 
0.1554,0.1558 and 0.1562 and give a ratio of pion to rho 
mass m^/rrip = 0.64, 0.59 and 0.50 respectively. We com- 
pare to unquenched results simulated using two flavors 
of Wilson fermions at re = 0.1575, re = 0.1580 |6j and 
k = 0.15825 .7] that give a ratio of pion to rho mass of 
0.69, 0.56 and 0.45 respectively. The nucleon isovector 
elastic form factors are extracted by solving the overcon- 
strained set of equations defined in Eq. (fTTTf) . 

In Fig. we compare results for the nucleon form fac- 
tors obtained in the quenched theory at re = 0.1554 using 
the local current and the lattice conserved current given 
in Eq. |5f . The local current is renormalized requiring 
charge conservation. The value of the renormalization 
constant that we find is G^ 1 (0) = Z v = 0.67. As we 
already pointed out, averaging over all directions of q, 
eliminates order a-terms in the lattice conserved current 
that arise from using fermionic fields at neighboring sites 
and the link variable that joins them. Despite the elim- 
ination of order a-terms, Fig. [S] shows a discrepancy be- 
tween results obtained using the local current and the 
lattice conserved current. In the case of the electric form 
factor, the results from the local current have smaller 
values, whereas for the magnetic case they have higher 
values. This means that these differences will be ampli- 
fied in the ratio of the two form factors. Since in our ap- 
proach the order a-terms are eliminated in the conserved 
current, there is no other obvious improvement, as far as 
the current is considered, that we can implement. Given 
that the lattice current is the one that ensures charge 
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FIG. 6: Ge and Gm as functions of Q 2 in the quenched theory 
at re = 0.1554 extracted using the local current ip(x)'y ll ij}(x) 
(open circles) and the symmetrized lattice conserved electro- 
magnetic current (crosses). 



conservation, the consistent approach is to use this cur- 
rent. A further argument for this choice is provided by 
considering the electric form factor Ge, which can be 
evaluated using Eq. (JSJ) or Eq. 10. In Fig. 0we compare 
the electric form factor extracted using Eq. © with that 
extracted from Eq. (JSJ. As can be seen for Q 2 larger 
than about 0.5 GeV 2 there is perfect agreement when 
using the lattice conserved current. The disagreement 
at lower Q 2 can be understood from the dependence on 
the momentum transfer appearing in the right hand side 
of Eq. ©. As q -> the right hand side of Eq. © 
tends to zero. Inverting to obtain Ge from the measured 
11(0, — q; T4; fx = i) becomes inaccurate resulting in an 
erroneous value for Ge- Using the local current on the 
other hand, we observed small differences between the 
results obtained using Eq. © and Eq. © up to values of 
Q 2 as large as 2 GeV 2 , which indicates a lattice artifact. 
Therefore, given charge conservation and consistency of 
the results coming from two different determinations of 
Ge, we conclude that the lattice conserved current is the 
best choice for the evaluation of the form factors within 
our current framework. Furthermore for Q 2 < 0.5 GeV 2 , 
we will only use Eq. for the determination of Ge 
whereas for higher Q 2 values both Eqs. |JSJ and @ will 
be used. 
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FIG. 7: Gb as a function of Q in the quenched theory at 
k — 0.1554 extracted using Eq. ||SJ (open circles) when the 
electromagnetic current is in the spatial direction and using 
Eq. © (crosses) with the electromagnetic current in the tem- 
poral direction. The top graph show results extracted using 
the local current tp(x)y f i,'il>(x) and the lower using the sym- 
metrised lattice conserved current. 



In Fig.|Hlwe show the results for the electric form fac- 
tor at three values of k for the quenched and the un- 
quenched cases. On the scale of this figure, only a weak 
quark mass dependence is seen. Both quenched and un- 
quenchcd results decrease as the quark mass decreases, 
yielding a larger slope at small Q 2 , which is the expected 
behavior. In order to better resolve differences in our 
data, we plot in Fig. [5] the ratio of the electric form 
factor to the proton dipole form factor, Gs/Gd- Both 
quenched and unqucnched results are clearly higher than 
the experimentally determined data, decreasing with the 
quark mass. The unquenched results in general show a 
stronger quark mass dependence leading to smaller val- 
ues in the chiral limit. The main observation, however, is 
that both quenched and unquenched results have a differ- 
ent Q 2 dependence as compared to the results extracted 
from experimental measurements: The lattice data have 
a positive slope at small Q 2 whereas experiment favors a 
negative slope. The two main uncertainties regarding the 
lattice results are finite a-effects and whether we are close 
enough to the chiral limit. Since unquenched Wilson con- 
figurations are only available at this lattice spacing, as- 
sessing whether finite a-effects can explain this behavior 
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FIG. 8: The isovector form factor Ge as a function of Q 2 . 
We show quenched lattice results at k = 0.1554 (crosses), at 
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unquenched results are shown at k — 0.1575 (open triangles) 
k — 0.1580 (filled circles) and at k, — 0.15825 (open squares). 
The filled triangles show experimental results for the isovector 
electric form factor extracted using the analysis described in 
Section III and data from Refs. HUE!- 



3.0 
2.5 
2.0 
1.5 
1.0 
0.5 
0.0 



X N F =0 m^O.56 GeV 
O N F =0 m"=0.49 GeV 
* N F =0 m"=0.41 GeV 
A experiment (isovector) 

A 

I „ 



• N„: 



=2 m =0.69 GeV 
=2 m =0.51 GeV 



NI=2 m™=0.38 GeV 

1 




0.0 



0.5 



1.0 

Q 2 (GeV 2 ) 



1.5 



2.0 



FIG. 9: The isovector electric form factor, Ge, divided by the 
proton dipole form factor as a function of Q 2 . The notation 
is the same as in Fig. |H| 



is beyond the scope of the present study. Also dynamical 
Wilson configurations at smaller quark masses on large 
enough volumes are not available so at present we cannot 
evaluate these form factors closer to the chiral limit. 

The evaluation of the magnetic form factor Gm is 
done using Eq. (|1U|> . which employs the optimal source. 
The results for the magnetic form factor are shown in 
Fig. ^| and, on the scale of this figure, the lattice re- 
sults are closer to experiment than the results for Ge- 
Again, quenched and unquenched results decrease with 
the quark mass with the unquenched results showing a 
stronger quark mass dependence. The stronger quark 
mass dependence of the unquenched data at low Q 2 is 
more clearly seen in Fig. 1111 where we plot the ratio 
Gm/G^. Again in the chiral limit we expected a reduc- 
tion in the value of Gm bringing lattice results closer to 
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FIG. 10: The isovector magnetic form factor, Gm, as a func- 
tion of Q 2 . The notation is the same as in Fig. |H1 
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FIG. 11: The isovector form factor, Gm, divided by the pro- 
ton dipole form factor taking fi — 4.71, as a function of Q 2 . 
The notation is the same as in Fig. |g| 



experiment. It is worth noting that the experimentally 
determined isovector form factor is very well described by 
the dipole form Gd(Q 2 ) whereas the lattice data clearly 
show deviations from the dipole form at least for the mass 
range considered in this work. To directly compare, how- 
ever, to experiment one has to carry out a chiral extrap- 
olation of Ge and Gm- This is discussed in the next sec- 
tion. Ideally one must also carry out the continuum limit 
using lattices of different values of a, which however is be- 
yond the scope of the present work. For completeness, 
we show in Fig.^lthe form factors Fx and F2, which are 
a linear combination of Ge and Gm ■ What can be seen 
is that, in the case of Fx, the lattice results show only a 
very weak increase in the slope as the quark decreases. 
The slope of Fx is directly related to the transverse size of 
the hadron |4J and one expects an increase in the slope 
as the quark mass decreases, which is not observed in 
the lattice data. On the other hand, in the case of F2 
one observes a stronger quark mass dependence. This 
stronger quark mass dependence potentially can lead to 
agreement with experiment after the chiral extrapolation 



1.0 



0.8 - 



0.6 



0.4 



0.2 



0.0 



2 - 



~i 1 1 

x N F =0 m^O.56 GeV 

o N F =0 m^O.49 GeV 

* N F =0 m n =0.41 GeV 




A N F =2 m n =0.69 GeV 
• N p =2 m n =0.51 GeV 

a N =2 m =0.38 GeV 

F TT 

A experiment (isovector) 



4. 



__i £ 



0.0 0.5 1.0 1.5 2.0 

Q 2 (GeV 2 ) 

FIG. 12: The isovector form factors Fx (top) and F2 (bottom), 
as a function of Q 2 . The notation is the same as in Fig. [H] 



2.0 



cj 1.5 
o 

1.0 



0.5 



0.0 



X N F =0 m^O.56 GeV 
O N F =0 m*=0.49 GeV 
* N F =0 m"=0.41 GeV 
A experiment (isovector) 



A N F =2 m^O.69 GeV 
• N =2 m =0.51 GeV 




0.0 0.5 1.0 1.5 2.0 

Q 2 (GeV 2 ) 



2.5 



FIG. 13: The isovector ratio, (j,Ge/Gm, as a function of Q 2 . 
The notation is the same as in Fig. [£] 



is carried out. In the case of Fx, given the weaker quark 
mass dependence and the larger deviation from experi- 
ment, one would require a non-trivial mass dependence 
at small quark masses to reconcile the lattice data with 
experiment. The experimentally interesting ratio of form 
factors, ijlGe/Gm, is shown in Fig. El As can be seen 
it shows very little dependence on the quark mass and, 
modulo finite a-effects, it can already be compared to ex- 
periment. The ratios obtained in the quenched and the 
unqucnchcd theory are in agreement with each other but 
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disagree with the behavior extracted from experiment. 
This disagreement is also clearly seen in the ratio F2/F1 
shown in Fig. 1141 especially at small Q 2 . 



V. EXTRAPOLATION TO THE CHIRAL LIMIT 

In order to compare our results for Ge and Gm with 
experiment, we must extrapolate the lattice results to the 
chiral limit. The quark masses used in this work corre- 
spond to pion masses in the range 560 to 410 MeV in the 
quenched theory and 690 to 380 MeV in the unqucnchcd 
theory. Pion cloud effects are expected to be small in this 
range of pion masses and therefore we expect the results 
to show a linear dependence in the pion mass squared, 
m 2 . To carry out correctly the chiral extrapolation of the 
form factors one would need chiral perturbation theory in 
the range of pion masses that we have results and valid for 
momentum transfers Q 2 up to about 2 GeV 2 . The only 
chiral expansion for the form factors presented recently is 
limited to small momentum transfers [2]]]. On the other 

hand, one expects that for values of Q 2 ~ 0.5 GeV 2 non- 
analytic terms are suppressed and a linear dependence in 
m 2 provides a good description to the data. In Figll5lwe 
plot the magnetic and electric form factors for the three 
lowest Q 2 values as a function of m 2 . We used Eq. © 
to extract the electric form factor since Q 2 < 0.5 GeV 2 
and only Eq. JUJ yields reliable results. For these low- 
est values of Q 2 non-analytic terms could be important 
and should be visible as the pion mass decreases. As 
can be seen in this figure, a linear dependence is consis- 
tent for the quenched data at the three lowest Q 2 values. 
The unquenched data are also consistent with a linear 
behavior. This is indeed what is observed also for the 
higher Q 2 values as can be seen in Fig. ^] where we show 
the unquenched results at Qq — 1.37 GeV 2 where by 
Qq = 2Mn(En — Mn) we denote the momentum transfer- 
squared in the chiral limit obtained by using the physical 
nucleon mass. It is therefore reasonable to extrapolate 



the form factors linearly in to 2 to obtain results in the 
chiral limit. Since Q 2 depends on the mass of the nucleon 
it changes with the quark mass and we need to extrapo- 
late form factors evaluated at somewhat different Q 2 . To 
leading order Q 2 decreases linearly with m 2 . Therefore 
we perform a fit to the form 



/(Q 2 ,m 2 ) = /(Q 2 ,0)+^m 2 



(19) 



where we extract the form factor at the chiral limit at Qq . 
As can be seen in Fig.ll5lthis linear behaviour is well sat- 
isfied for Qg = 1.37 GeV 2 the largest value shown in the 
figure. Another option is to interpolate lattice data at 
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FIG. 15: The isovector electric (upper) and magnetic(lower) 
form factor as a function of m 2 for the three lowest Q 2 values 
available on our quenched lattice. With Q% we denote the 
momentum transfer square in the chiral limit. Unquenched 
results are shown with the filled circles at the lowest available 
Q 2 -value and with the open triangles at Q 2 , = 1.37 GeV 2 . 

Qg and perform a chiral extrapolation at this constant 
value of the momentum transfer squared. Only the value 
A of the slope should be affected whereas f(Qo, 0) should 
not change. We checked these two procedures for vari- 
ous values of Q 2 . We find that the results for /(QqjO) 
obtained using these two procedures are indeed consis- 
tent. In what follows we will therefore use Eq. i|19|) for 
the chiral extrapolation. 

The resulting values at the chiral limit are shown in 
Fig. El for Ge and Gm- The disagreement with experi- 
ment is larger in the case of the electric form factor and 
can be traced to its weak quark mass dependence. Lat- 
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tice results also show a different Q 2 dependence as com- 
pared to experiment. The linearly extrapolated lattice 
results are closer to experiment in the case of the mag- 
netic form factor. For comparison, we show in Fig. 1171 
the lattice results for F\ and F2 after linearly extrapo- 
lated to the chiral limit. There is little deviation between 
unquenched and quenched results at the chiral limit. In 
addition, for Q 2 > 0.5 GeV 2 there is good agreement be- 
tween the lattice results for F2 and the results extracted 
from experiment. This is not the case for F% where the 
experimentally determined isovector F\ decays faster as 
compared to the lattice results. 
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FIG. 17: The form factors Fi (upper) and F2 (lower) as a 
function of Q 2 . Crosses show quenched results at the chi- 
ral limit, and filled circles unquenched results. The results 
extracted from experiment are shown by the filled triangles. 



G E (Q 2 ) = 



1 



1+01 



(21) 



0.0 ' ' 1 1 1 

0.0 0.5 1.0 1.5 2.0 

Q 2 (GeV 2 ) 

FIG. 16: The ratios Ge/G c i (upper) and ^lGm/Gci (lower) 
as a function of Q 2 at the chiral limit. Open triangles show 
quenched results obtained at the chiral limit by linear ex- 
trapolation of the form factors Gm and Ge and filled circles 
denote the corresponding unquenched results. The results for 
these isovector ratios extracted from experiment are shown 
by the filled triangles. 

In order to obtain the isovector magnetic moment 
Gm(0), one needs to evaluate the magnetic form factor 
Gm at Q 2 = 0. This requires an extrapolation of lattice 
results to Q 2 = 0. One fitting Ansatz commonly used 
to describe the Q 2 dependence of the form factors is a 
dipole. We thus consider a dipole form with different 
isovector magnetic and electric dipole masses squared, 
M m and M e : 

G M (Q 2 ) = / Gm( ° ) n2 (20) 



A good fit giving x 2 / (degree of freedom) ~ 1 is obtained 
when the form factors are fitted for Q 2 ~ 2.5 GeV 2 . 
The quality of the fits is shown in Fig. 1181 where we 
show quenched and unquenched data at n = 0.1558 and 
K = 0.1580 respectively. On the same figure, we include 
a fit to an exponential form. For the magnetic form fac- 
tors, an exponential Ansatz describes very well the Q 2 
dependence, especially in the quenched case where it is 
in fact favored. An exponential fit, however, does not 
provide a good fit to the electric form factors. Therefore 
in order to extract the r.m.s. radii we will use through- 
out a dipole Ansatz. In Table [H] we give the magnetic 
moment, Gm(0), and the magnetic and electric dipole 
masses extracted from the dipole fits to the quenched 
and unquenched results at each value of the quark mass. 
The dipole masses extracted from the fits are generally 
larger than the value of Md = 0.71 GeV 2 entering in the 
proton form factor Gd- This is consistent with the fact 
that the lattice data normalized with Gd are not constant 
but increase as a function of Q 2 as shown in Figs. [H] and 
im In the same Table we also give the values extracted 
by applying dipole fits to the form factors after they have 
been linearly extrapolated to the chiral limit. Since the 
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dipole masses decrease with the quark mass, the fits of for M e and M m . In the unqucnchcd case these values are 
the form factors at the chiral limit yield a smaller value rather close to the value of the proton dipole mass, M<j. 



TABLE II: 





Gm(0) 


M m (GeV 2 ) 


M e (GeV 2 ) 


< r? > 1/2 (fm) 


<r 2 2 > 1 / 2 (fm) 


0.1554 
0.1558 
0.1562 
k c =0.1571 


Quer 
4.11(7) 
4.02(8) 
3.90(9) 
3.73(13) 


iched 32 3 x ( 
1.29(4) 
1.28(4) 
1.19(4) 
1.03(5) 


34 a -1 
1.24(1) 
1.15(1) 
1.08(1) 
0.90(2) 


= 2.14(6) GeV 
0.520(5) 
0.538(6) 
0.550(8) 
0.585(13) 


0.64(1) 
0.64(1) 
0.66(1) 
0.72(2) 


0.1575 
0.1580 

0.15825 
Kc =0.1585 


Unqu< 
4.45(14) 
4.34(43) 

Unquc 
4.10(46) 
3.25(48) 


;nched 24 3 x 
1.53(7) 
1.23(16) 

inched 24 3 x 
1.17(17) 
0.792(17) 


40 a' 1 
1.55(1) 
1.41(2) 

32 a' 1 
1.19(4) 
0.66(4) 


= 2.56(10) GeV 
0.467(7) 
0.462(23) 

= 2.56(10) GeV 
0.500(29) 
0.756(36) 


0.58(2) 
0.67(5) 

0.68(6) 
0.79(13) 



The quark mass dependence of the magnetic moments 
extracted from the dipole fits is presented in Fig. A 
linear extrapolation of the quenched results to the chiral 
limit leads to Gm(0) = 3.67(3). As expected, this value 
is in agreement with the value obtained from the dipole 
fit of Gm at the chiral limit quoted in Table ITT1 This 
agreement confirms that our extrapolation using Eq. I|19|) 
yields results consistent to those obtained by extrapo- 
lating to the chiral limit at constant value of Q 2 . The 
slope of Fi at Q 2 = determines the transverse size of 
the hadron, < r\ >— —AdFi/ dQ 2 \q2 =() . In the non- 
relativistic limit the r.m.s. radius is related to the slope 
of the form factor at zero momentum transfer. Therefore 
the r.m.s. radii can be directly obtained from the values 
of the dipole masses by using 



< 



rf>- 



6 dF^Q 2 ). 



12 

M 



= ttt i= 1,2 



Fi{Q 2 ) dQ 2 ,w 

(22) 

The electric and magnetic radii are given by < r 2 m >= 
12/M e m and can be directly evaluated from the values 
given in Table [H] We can also obtain < rf > in terms of 
M m and M e using the relations 



TABLE III: The first column lists the fixed parameters and 
the second their values at the physical pion mass. The third 
column gives the fitted parameters and the fourth the values 
determined by fitting to Eq. 1241 for the magnetic moment, 
to Eq. for r\ and to Eq. J27| for r\. 



3F 2 (0) 



< r, >= 



12(1 + F 2 (0)) 



2 _ 12 

<ri> "Me 2M 2 N 1 F 2 (Q)M m F 2 (0) 

(23) 

In Eq. |g3J| we take F 2 (0) = G M (0) - 1 extracted from 
the dipole fits. Alternatively, one can fit directly the F-y 
and F2 form factors and obtain the dipole masses Mi 
and M2. What one finds via this procedure is that the 
values for < rf > tend to be smaller than but consis- 
tent within errors with the ones extracted using Eq. (|23|l . 
The quark mass dependence of the magnetic and elec- 



Fixed 


Empirical value 


Fitted 


Fitted value 


parameter 




parameter 




9A 


1.267 


Kv(0) 


6.08(11) 


CA 


1.125 


cv 


-2.75(50) GeV" 1 


F-jt 


0.0924 GeV 


E! 


-5.60 (5) GeV -3 


M N 


0.9389 GeV 


B w 


-0.3(3) GeV -3 


A 


0.2711 GeV 


B C 2 


0.61(4) 



trie dipole masses is also shown in Fig. A linear 
dependence in m 2 is consistent for the quenched results 
yielding, at the chiral limit, M m = 1.09(10) GeV 2 and 
M e = 0.89(4) GeV 2 . Again, these values are in agree- 
ment with the values extracted by fitting the form factors 
after they have been linearly extrapolated to the chiral 
limit as can be seen from the values quoted in Table ITU 
After scaling the magnetic moments by the ratio of the 
physical nucleon mass to the one measured on the lat- 
tice, they become an increasing faction of m 2 as can be 
seen in Fig. ^5] In the same figure we also include the 
quenched results obtained from Ref. [5j, that used per- 
turbatively improved Wilson fermions. We choose data 
at j3 = 6.0 where the value of the lattice spacing ex- 
tracted from their nucleon mass is a -1 = 1.83 GeV. This 
is close enough to our quenched lattice to allow a mean- 
ingful comparison. It is reassuring that, despite the fact 
that different currents and Wilson fermions were used in 
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FIG. 18: The magnetic (upper) and electric (lower) form fac- 
tors for the quenched and the unquenched cases with fits to 
dipole (dashed line) and exponential (solid line) forms. 



the two calculations, results at a similar pion mass are 
consistent. 

As the pion mass decreases, one expects cloud pion 
contributions to become important and deviations from 
the linear dependence on should be observed, in par- 
ticular at low Q 2 , thus affecting the values of Gm(0) and 
the dipole masses and hence the r.m.s. radii. In a recent 
calculation the quark mass dependence of the isovector 
magnetic moment and radii was determined. This was 
done within a chiral effective theory with explicit nucleon 
and A degrees of freedom yj, |22( . The isovector anoma- 
lous magnetic moment as a function of the pion mass to 
one-loop order is given by [2^ 
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FIG. 19: The magnetic moment (upper), the magnetic (mid- 
dle) and electric (lower) dipole mass extracted from fits, as- 
suming a dipole dependence for the form factors, are shown 
as a function of the pion mass squared. The dashed line is 
a linear fit to the quenched results. The values at the chiral 
limit are obtained from fitting the linearly extrapolated form 
factors at the chiral limit to a dipole Ansatz. Our quenched 
data are shown by the crosses, the unquenched data by the 
filled circles, our scaled quenched (unquenched) magnetic mo- 
ments by the open triangles (open triangles) and the quenched 
results at /3 = 6.0 of Ref. _5] by the asterisks. 




where 



Ri(m) 



VA 2 



2A 



log 



A + ^fA 2 



A-VA 2 ^ 



(24) 



(25) 

and A = Ma — Mn is the A-nucleon mass splitting. Fol- 
lowing Ref. [22| we fix g A , c A) F„, Mn and A to their 
physical values given in Table ITTT1 and vary k v (0), cy and 
Ei . The counter term E\ depends on the regularization 
scale A for which we take A = 0.6 GeV in order to make 
contact with Ref. j2^|. As can be seen in Fig.^|the mag- 
netic moments in the quenched and unquenched theory 
are in agreement and therefore we use both sets to fit to 
the chiral effective theory result given in Eq. (|24|l . Fit- 
ting to the rescaled data we obtain the curve shown by 
the solid line in Fig.[3U| In Table ITTT1 we give the values 
of k v (0), cy and E\ extracted from the fit. The dashed 
lines give the maximal error band determined by varying 
the fitted parameters by the quoted errors. The extrapo- 
lated value of the magnetic moment at the physical pion 
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FIG. 20: Chiral extrapolation of the magnetic moment (up- 
per) and the r.m.s radii r\ (middle) and r% (lower). The solid 
line is the best fit to the effective chiral theory results. The 
dashed lines show the maximal allowed error band using the 
errors on the fitted parameters. 

mass is in agreement with experiment. 

For the isovector Dirac and Pauli radii we use the one- 
loop results given in Ref. [j|: 
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12Bio 



1 
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,-2 



(lOg 2 A + 2) log (^) 



(47TfV) 2 54tt 2 F2 
+30R 2 {m v ) 



26 + 30 log (— ^ 



(26) 



and 



1 



( g\M N 



K v (mi) [ 8F 2 7rm T 
+2AM N B c2 



c\M N 



(27) 



where 



Ba(m) 



2VA 2 



= log 



A + VA 2 



A- VA2 



(28) 

The only parameter that we vary in fitting the Dirac ra- 
dius is the counter term Bio, which depends on the scale 



A and parametrizes short distance contributions. Once 
the magnetic moment is fitted, the only parameter en- 
tering in the Pauli radius that we vary is the counter 
term B C 2, which is the analogue of £>io. The resulting 
fits for the radii are shown in Fig. [2UJ The pion mass 
dependence of the Dirac radius is not well reproduced. 
Since this is related to the slope of F% this is not surpris- 
ing given that the lattice results have a different slope 
from the experimental one and hardly show any quark 
mass dependence. 



VI. CONCLUSIONS 

The elastic isovector nucleon form factors are calcu- 
lated in lattice QCD with Wilson fermions both in the 
quenched approximation and using unquenched config- 
urations [|| 0] with two flavors of dynamical Wilson 
fermions. The current work presents an improvement to 
a previous lattice study |5| , carried out in the quenched 
approximation, in a number of ways: In the quenched 
theory we use a lattice of twice the spatial and tempo- 
ral size. This allows an accurate determination of the 
form factors at lower values of Q 2 , enabling us to ex- 
tract more reliably the dependence on Q 2 . In addition, 
the quenched calculation is carried out at smaller quark 
masses, bringing us closer to the chiral limit. Preliminary 
quenched results on a lattice of size 32 3 x 48 at = 6.0 
using Wilson fermions were presented in Ref. |23| . Al- 
though the low Q 2 range probed is the same as in the 
current work, only results after linear extrapolation to 
the chiral limit, using data computed at two light quark 
masses and thus carrying large statistical errors, were 
discussed. Furthermore in this work we evaluated the 
form factors in the unquenched theory allowing us to as- 
sess unquenching effects for pion masses down to about 
380 MeV. Finally, an improved overconstrained analy- 
sis is carried out where the nucleon source is optimized 
and all the lattice momentum vectors contributing to a 
given value of Q 2 are taken into account. The result- 
ing statistical errors are therefore small enough that a 
comparison between quenched and unquenched results is 
meaningful. What we find is that both quenched and 
unquenched results for both form factors decrease with 
the quark mass. Unquenching effects are small and the 
results obtained after a linear extrapolation in m 2 to the 
chiral limit fall on the same curve as can be seen for exam- 
ple in Fig. E|f° r Fi an d F%. Assuming lattice artifacts 
are under control, this improved analysis gives results 
that can be compared to experiment by extracting the 
isovector form factors from the proton and neutron mea- 
surements of these quantities. The largest uncertainties 
regarding our lattice results are how close to the contin- 
uum limit these results are and the chiral extrapolation. 
With these caveats in mind, the comparison of the results 
obtained here to experiment reveals interesting features: 
Both quenched and unquenched results are higher than 
the experimentally extracted form factors, with the de- 
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viations being larger in the case of the electric isovector 
form factor. In the quenched case, where we have very ac- 
curate results at low momentum transfer, we find that the 
electric form factor decreases slower with Q 2 compared to 
what is observed experimentally. This different behavior 
is also reflected in the ratio ^.Ge/Gm, where the lattice 
results are constant up to about Q 2 = 2.5 GeV 2 whereas 
the experimentally determined data decrease as a func- 
tion of Q 2 . This leads to smaller values for the r.m.s. 
radii showing that pion cloud contributions are smaller 
at these quark masses than in the physical regime. In 
the range of quark masses investigated in this work, the 
quark mass dependence observed for [iGe/Gm is small 
and so are unqucnching effects. Using chiral effective 
theory to one-loop to extrapolate the magnetic moment 
to the chiral limit, we find that the lattice results ex- 
trapolate nicely to the experimental value. The charge 
radius, on the other hand, is constant over the range of 
quark masses used in this work and therefore deviates 
from experiment. This is again related to the deviation 
observed between experiment and lattice results in the 
case of the charge form factor. Since the lattice size in 
the quenched case is large, we do not expect finite vol- 
ume effects to be the reason for the discrepancy. What 
needs to be checked is finite lattice spacing effects and 
whether we are close enough to the chiral limit. The ob- 
served disagreement with experiment in the case of the 
charge form factor is puzzling and a study using finer 
lattices should follow. If one requires in addition dynam- 



ical fermions and small quark masses to be closer to the 
chiral limit, then such a study would require large com- 
puter resources. Unquenched configurations with pion 
masses down to 250 MeV on reasonably large and fine 
lattices will become available in the near future enabling 
us, using the techniques of the current work, to obtain re- 
sults closer to the physical regime avoiding uncontrolled 
extrapolations. 

Acknowledgments: We would like to thank B. Orth, 
Th. Lippert and K. Schilling [6( as well as C. Urbach, 
K. Jansen, A. Shindler and U. Wenger for providing 
the unquenched configurations used in this work. A.T. 
would like to acknowledge support by the University of 
Cyprus and the program "Pythagoras" of the Greek Min- 
istry of Education. C.A. and J.W.N, would like to thank 
the Institute for Nuclear Physics at the University of 
Washington for its hospitality and partial support dur- 
ing the completion of this work. The computations for 
this work were partly carried out on the IBM machine at 
NIC, Julich, Germany. This work is supported in part 
by the EU Integrated Infrastructure Initiative Hadron 
Physics (I3HP) under contract PJI3-CT-2004-506078 and 
by the U.S. Department of Energy (D.O.E.) Office of Nu- 
clear Physics under contract DF-FC02-94ER40818 This 
research used resources of the National Energy Research 
Scientific Computing Center, which is supported by the 
Office of Science of the U.S. Department of Energy under 
Contract No. DE-AC03-76SF00098. 



[1] M. Jones, et al. Phys. Rev. Lett. 84, 1398 (2000); O. 

Gayou, et al, Phys. Rev. C 64, 038202 (2001); O. Gayou 

et al, Phys. Rev. Lett. 88, 092301 (2002). 
[2] H.-Y. Gao, Int. J. Mod. Phys. E 12, 1 (2003). 
[3] C. E. Hyde-Wright and K. de Jager, Ann. Rev. Nucl. 

Part. Sci. 54, (2004) 217. 
[4] M. Burkardt, Int. J. Mod. phys. A 18, 173 (2003). 
[5] M. Gockeler et. al, Phys. Rev. D 71, 034508 (2005). 
[6] B. Orth et. al, Phys. Rev. D 72, 014503 (2005). 
[7] C. Urbach et. al, Comput. Phys. Commun.174, (2006) 

87. 

[8] D.B. Leinweber and R. M. Woloshyn, Phys. Rev. D 43, 
(1991) 1659. 

[9] LHP, SESAM Collaboration, Phys. Rev. D 66, 034506 
(2002). 

[10] C. Alexandrou et al. Phys. Rev. Lett. 94, 021601 (2005). 
[11] C. Alexandrou et al. Phys. Rev. D 69, 114506 (2004). 
[12] C. Alexandrou et. al, Nucl. Phys. B414, (1994) 815. 
[13] M. Albanese et. al, Phys. Lett. B 192, 163 (1987). 
[14] C. Morningstar and M. J. Peardon, Phys. Rev. D 69, 
054501 (2004). 

[15] A. Hasenfratz and F. Knechtli, Phys. Rev. D 64, 034504 
(2001). 

[16] L. Andivahis et al, Phys. Rev. D 50 5491 (1994); R. C. 



Walker et al, Phys. Rev. D 49 5671 (1994); A. F. Sill, 
et al, Phys. Rev. D 48 29 (1993); P. E. Bosted, et al, 
Phys. Rev. C 42 38 (1990); F. Borkowski, et al, Nucl. 
Phys. B93, 461 (1975). 
[17] H. Anklin, et al, Phys. Lett. B 336, 313 (1994); H. 
Anklin, et al, Phys. Lett. B 428, 248 (1998); G. Kubon, 
et al, Phys. Lett. B 524, 26 (2002); W. Xu, al. Phys. 
Rev. Lett. 85 2900 (2000); W. Xu, et al, Phys. Rev. C 
67, 012201 (2003); A. Lung, Phys. Rev. Lett. 70, 718 
(1993); R. Madey, et al, Phys. Rev. Lett. 91, 122002 

(2003) ; G. Warren, et al, Phys. Rev. Lett. 92, 042301 

(2004) ; D. Rohe, et al, Phys. Rev. Lett. 83, 4257 (1999); 
J. Bermuth, et al, Phys. Lett. B 564 199 (2003). 

[18] S. Galster et. al, Nucl. Phys. B32,(1971) 221. 

[19] B. D. Milbrath, et al. Phys Rev. Lett. 80 452 (1998); 

Th. Pospischil, et al. Eur. Phys. J. A12 125 (2001); V. 

Punjabi et al. Phys. Rev. C 71 055202 (2005). 
[20] S. Necco and R. Sommer,Nucl. Phys. B622 328 (2002). 
[21] S. Ando, J.-W. Chen and C.-W. Kao, hep-ph/0602200 
[22] T. R. Hemmert and W. Weise, Eur. Phys. J. A 15,487 

(2002). 

[23] T. Bakeyer et al, (QCDSF-UKQCD Collaboration), 
Nucl. Phys. Proc. Suppl.128, 82(2004). 



